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ABSTRACT 

We re-address the classification criterion for extremely red galaxies (ERGs) of Pozzetti 
& Mannucci (2000 - PM00), which aims to separate, in the l c — K (or R c — K) vs. 
J — K colour-colour diagram, passively evolving, old (^ 1 Gyr) stellar populations in 
a dust-free environment, associated with ellipticals (Es), from dusty starburst galaxies 
(DSGs), both at 1 < z < 2. We explore a category of objects not considered previously, 
i.e., galaxies forming in this redshift range on short (0.1 Gyr) timescales and observed 
also in their early, dusty post-starburst phase. We also investigate the impact of struc- 
ture of the dusty medium and dust amount on the observed optical/near-IR colours of 
high-z DSGs/DPSGs, through multiple-scattering radiative transfer calculations for 
a dust/stars configuration and an extinction function calibrated with nearby dusty 
starbursts. As a main result, we find that dusty post-starburst galaxies (DPSGs), 
with ages between 0.2 and 1 Gyr, at 1.3 < z < 2 mix with Es at 1 < z < 2 for a 
large range in dust amount. This "intrusion" is a source of concern for the present 
two-colour classification of ERGs. On the other hand, we confirm, in agreement with 
PM00, that DSGs are well separated from Es, both at 1 < z < 2, in the I c — K vs. J — K 
colour-colour diagram, whatever the structure (two-phase clumpy or homogeneous) 
of their dusty medium and their dust amount are. This result holds under the new 
hypothesis of high-z Es being as dusty as nearby ones. Thus the interpretation of the 
optical/near-IR colours of high-z Es may suffer from a multiple degeneracy among 
age, metallicity, dust and redshift. We also find that DPSGs at z <~ 1 mix with DSGs 
at 1 < z < 2, as a function of dust amount and structure of the dusty medium. All 
these results help explaining the complexity of the ERG classification emerged from 
recent surveys. Model DSGs/DPSGs with K < 20.5 and selected as ERGs on the basis 
of their observed optical/near-IR colours span a large range in absorption efficiency 
and/or intrinsic UV-to-optical luminosity ratio. As a consequence, the total emission 
from dust in the rest-frame mid-IR-(sub)mm ranges from ~ 10 10 to ~ 5.6 x 10 12 Lq 
for DPSGs through DSGs. This large dynamic range is consistent with estimates made 
from submm observations of bright (K < 19.5-20) ERGs. Finally, we show that ERG 
samples are biased towards a lower fraction either of DSGs/DPSGs at 1.3 < z < 2 or 
of DSGs/DPSGs at z ~ 1, when selection is made via either the R c — K 5.3 or the 
I c — K ^ 4 colour criterion. 

Key words: radiative transfer - dust, extinction - galaxies: elliptical and lenticular, 
cD - galaxies: starburst - galaxies: high redshifts. 



1 INTRODUCTION 

The extremely red objects (EROs), originally discovered in 
the K-band survey of Elston, Rieke & Rieke (1988), have 
become more and more present in deep near-IR and optical 
surveys (see Yan & Thompson 2003 for a review). Oper- 
ationally defined in terms of the observed optical/near-IR 
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colours (e.g. R c — K ^ 5.3 or I c — K ^ 4), significantly differ- 
ent from those of typical field objects, EROs may neverthe- 
less comprise an assorted mix of Galactic and extragalactic 
objects. We focus on the extremely red galaxies (ERGs), 
which account for the bulk of the observed EROs. For the 
optically selected ERGs, the contamination from AGNs is 
less than 15 per cent (Alexander et al. 2002; Brusa et al. 
2002; Smail et al. 2002), even though obscured, Compton- 
thin AGNs, partly associated with dusty starbursts having 
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large star formation rates (SFRs), most probably dominate 
X-ray selected ERG samples (Stevens et al. 2003). 

ERGs still lack an uncontroversial classification because 
their very red colours can be reproduced by at least two 
very different sets of spectral energy distributions (SEDs) 
not associated with an AGN: those expected for passively 
evolving, old 1 Gyr) stellar populations in a dust-free 
environment, associated with ellipticals, and those expected 
for dusty starburst galaxies (DSGs), both at high redshift 
(z > 1). Indeed it emerges that a half or more of the entire 
ERG sample in the literature contains ellipticals (Moriondo, 
Cimatti & Daddi 2000; Stiavelli & Treu 2001; Mannucci et 
al. 2002; Cimatti et al. 2002a). 

ERGs represent a key-case in the knowledge of the cos- 
mic history of star formation. In fact old, massive ellipticals 
at z > 1 represent a problem for the standard hierarchical 
merging scenario of galaxy formation and evolution (White 
& Frenk 1991), where the most massive galaxies are assem- 
bled through a merging tree of less massive systems, even 
though the bulk of their stars may already be very old at 
the time of the last major merging (e.g. Kauffmann, Chariot 
& White 1996; Springel & Hernquist 2003). Conversely, the 
association of ERGs with DSGs at intermediate redshifts 
does not necessarily favour the latter scenario, if galaxies 
are assembled within times comparable to their star for- 
mation timescales and the latter are shorter the larger the 
masses of the galaxies are, as independently proposed for 
late- (Gavazzi, Pierini & Boselli 1996) and early-type galax- 
ies (Thomas, Maraston & Bender 2002). Interestingly "hy- 
per" ERGs have been interpreted as primodial ellipticals in 
a dusty starburst phase at z ~ 3 (Totani et al. 2001). 

Unfortunately the classification either as a passively 
evolving, old elliptical or as a dusty starburst galaxy at 
z > 1 is uncomfortably ambiguous in absence of an ex- 
tensively observed SED (e.g., Smith et al. 2001), but the 
needed observations are extremely challenging at present, 
given the faintness of the ERGs (e.g. Cimatti et al. 1999). 
Hence Pozzetti &: Mannucci (2000 - hereafter referred to 
as PM00) have proposed a classification criterion based on 
the relatively easily obtainable optical/near-IR photometry, 
particularly suitable for large samples of EROs with K < 20- 
21. This classification criterion is based on the rather neat 
difference in the distribution of the PM00 model ellipticals 
and DSGs at 1 < z < 2 in the I c - K (or R c - K) vs. 
J — K colour-colour diagram. This separation is mainly due 
to the fact that passively evolving, old stellar populations 
at these redshifts have their strong 4000 A break redshifted 
to A < 1.2 fim (J band), while coeval dusty galaxies show 
smoother SEDs and therefore redder observed J — K colours, 
whatever the dust/stars configuration is. 

For the dusty ERGs, there is no unanimous consensus 
on the typical dustiness and star formation rate. The detec- 
tion of ERO J164502+4626.4 (z = 1.44) at 850 ^m (Cimatti 
et al. 1998) suggested that dusty ERGs may host massive 
quantities of dust accompanied by very high SFRs, as in 
ultraluminous infrared galaxies (ULIRGs) (see also Cimatti 
ct al. 1997; Elbaz et al. 2002). Independently, 1.4 GHz ra- 
dio maps reveal that > 16 ± 5 per cent of the ERGs with 
K < 20.5 have radio luminosities expected for ULIRGs at 
z ~ 1 (Smail et al. 2002). Conversely, Andreani et al. (1999) 
and Mohan et al. (2002) conclude that about 20 per cent of 
the bright (K < 19.5-20) ERGs show strong submm emis- 



sion but that the dominant population of dusty, bright ERGs 
is probably made of starburst galaxies with a total mid-IR- 
mm rest-frame luminosity lower than 10 12 L0 and SFR < 
few xlO 2 M© yr" 1 (i.e., non-ULIRGs). 

Here we re-address the optical/near-IR photometric 
classification of ERGs and investigate the relation between 
reddening, dustiness and SFR in dusty ERGs. Four are the 
main new points of this study. First, we consider the pres- 
ence of stellar populations which are formed in a burst of 
short (0.1 Gyr) duration and are observed at 1 < z < 2 in 
an early post- starburst phase while still embedded in a dusty 
medium. Second, we consider the presence of passively evolv- 
ing, old, dusty elliptical galaxies at 1 < 2 < 2 as additional 
candidate ERGs. Third, the treatment of dust attenuation 
follows the Monte Carlo simulations of the radiative transfer 
by Witt & Gordon (2000 - hereafter referred to as WG00) 
for starbursts and post-starbursts, and those of Witt, Thron- 
son & Capuano (1992 - hereafter referred to as WTC92) 
for ellipticals. For DSGs and dusty post-starburst galaxies 
(DPSGs), these models explore the dependence of dust at- 
tenuation from wavelength as a function of a dust/stars con- 
figuration, structure of the dusty medium and dust proper- 
ties which were not considered simultaneously or at all by 
PM00 (see Sect. 2.2.2). Fourth, we calculate the total emis- 
sion from dust in the mid-IR-(sub)mm wavelength range for 
DSGs/DPSGs. 



2 MODELLING PASSIVELY EVOLVING, OLD 
ELLIPTICAL GALAXIES AND DUSTY 
STARBURST /POST-STARBURST 
GALAXIES AT 1< Z < 2 

For our candidate ERGs (passively evolving, old ellipticals 
and DSGs/DPSGs), we consider three representative red- 
shifts (z = 1.07, 1.39 and 1.85), consistent with published 
results. We adopt the so-called "concordance" cosmologi- 
cal model (fi m = 0.3, fl\ = 0.7) with Hubble constant 
Ho — 70 km s _1 Mpc -1 to link redshift and age of the 
models. 

We use two evolutionary synthesis codes (one for ellip- 
ticals and one for DSGs/DPSGs) and fix metallicity (i.e., 
solar - Zq) and initial mass function (Salpeter 1955), since 
PM00 made an exhaustive investigation of the influence of 
these parameters on the colours of their candidate ERGs. 
However we will discuss the impact of sub-solar metallici- 
ties on our results. 

Our focus is on the I c — K, R c — K and J — K colours, 
where Cousins (1978) I c and R c filters and Bessel & Brett 
(1988) J and K filters are used. 

2.1 Passively evolving, old elliptical galaxies 

2.1.1 Stellar populations 

As in PM00, passively evolving, old (^ 1 Gyr) elliptical 
galaxies are simulated through simple stellar population 
(SSP) models. Hereafter they are referred to as Es for sim- 
plicity. Their intrinsic (i.e., unattenuated by dust) SEDs are 
computed with the evolutionary synthesis code of Maraston 
(1998). 

We assume a maximal formation redshift Zf = 10, 
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against the three values of formation redshift (i.e., Zf = 3, 
4 and 6) assumed by PMOO for their model Es. An elliptical 
galaxy formed at 2/ = 10 has an age of about 5, 4 and 3 
Gyr at 2 = 1.07, 1.39 and 1.85, respectively. One formed 
at Zf = 3 (PMOO; see also Totani et al. 2001) has an age 
of about 3.4, 2.4 and 1.4 Gyr at 2 = 1.07, 1.39 and 1.85, 
respectively. 

We note that a 1 Gyr-old galaxy at z = 1.07, 1.39 or 
1.85 was formed at 2/ = 1.39, 1.85 or 2.59, respectively. Thus 
a 1 Gyr-old elliptical at 1 < z ^ 1.5 may have its progenitor 
in the redshift range that we are interested in. Hence it is 
important to investigate the case of a post-starburst galaxy 
and illustrate the difference between a model DPSG and 
a model elliptical galaxy of similar young/intermediate age 
(see Sect. 2.2.1). 

For model Es we adopt a minimum stellar mass of 7.5 x 

10 10 Mq, which is about the residual (i.e., after 25 per cent 
mass-loss owing to stellar evolution) stellar mass of a 1 Gyr- 
old SSP (from Maraston 1998) as well as of a 1 Gyr-old 
model DSG/DPSG (cf. Sect. 2.2.1), both originating from a 

10 11 Mq cloud of primordial gas. 

2.1.2 Dust attenuation 

We test the impact of dust attenuation on the observed 
SEDs of Es at 1 < z < 2. The residual interstellar medium of 
a dusty elliptical is considered to be negligible with respect 
to its characteristic stellar mass. 

For dusty Es we adopt the dust attenuation curves ob- 
tained from the Monte Carlo simulations of the radiative 
transfer by WTC92 for the CLOUDY geometry (see also 
WG00). This geometry assumes a spherical stellar distribu- 
tion that declines from the core as r and a homogeneous 
Milky Way-type dust present only in the inner portions of 
the model, extending from the center to two thirds of the 
galaxy radius and following a distribution that declines as 
r _1 . The results of this model are robust when different ra- 
dial dependences of the stellar and dust distributions are 
assumed (cf. Bianchi, Ferrara & Giovanardi 1996; Wise & 
Silva 1996; WG00). 

The value of the radial extinction optical depth from the 
center to the edge of the dust environment at V-band, as- 
suming a constant density, homogeneous distribution, (tv) 
is set to 0.5, consistent with the mean ry estimated from 
the analysis of broad-band colour gradients in nearby Es 
(WTC92; Wise & Silva 1996). These studies provide upper 
limits to Tv, since they do not take into account intrinsic 
gradients in metallicity of the stellar populations of Es (e.g. 
Saglia et al. 2000). 

In Fig. 2, we reproduce the attenuation function 1 , nor- 
malized to its V-band value, (A x norm ) of dusty Es for 
tv = 0.5, which produces a value of the stellar reddening 
(E(B - V)star) equal to 0.15. 



We stress the difference between dust extinction curve and dust 
attenuation curve. The former describes the absorption and scat- 
tering properties of a mix of dust grains of given size distribution 
and chemical composition as a function of wavelength; the latter 
is the convolution of the extinction curve with the geometry of the 
dusty stellar system. The attenuation curve depends on structure 
of the dusty medium and dust amount. 



2.2 Dusty starburst/post-starburst galaxies 

2.2.1 Stellar populations 

Intrinsic SEDs of young 0.1 Gyr)/intermediate-age (0.1- 
1 Gyr) starburst and intermediate-age post-starburst galax- 
ies are computed with the evolutionary synthesis code (ver- 
sion 2.0) PEGASE (Fioc & Rocca-Volmerange 1997), be- 
cause this code includes nebular emission. PEGASE 2.0 does 
not include the TP-AGB phase that is dominant at interme- 
diate ages, but this fact does not represent a major problem 
since we are considering observed wavelengths shorter than 
2.2 fim (K band). 

Gas is assumed to be transformed into Zq stars at a 
constant rate over a period tturst of 0.1 Gyr {short burst) or 
1 Gyr (long burst) at 1 < 2 < 2. For a total mass of trans- 
formed gas M gas = 10 11 Mq, a SFR equal to M gas /tburst 
corresponds to 10 2 Mq yr _1 (long burst) or 10 3 Mq yr _1 
(short burst). Note that SFRs ~ 10 3 Mq yr _1 have been 
inferred from observations of high-2 SCUBA submm sources 
(e.g., Smail, Ivison & Blain 1997; Hughes et al. 1998; Barger 
et al. 1998; Cimatti et al. 1998; Eales et al. 1999; Bertoldi 
et al. 2000; Lutz et al. 2001; Smail et al. 2003). 

These star formation histories and SFRs are comple- 
mentary to those adopted by PMOO for their DSGs for the 
following reason. The constant star formation rate from a 
formation redshift 2.5 or, especially, 6, assumed by PMOO 
for their basic model of a DSG, better corresponds to a con- 
tinuously star-forming galaxy with an intermediate/low SFR 
at 1 < 2 < 2. An object formed at 2/ = 2.5(6) is about 
0.7(2.3) Gyr old at 2 = 2, and 3.2(4.8) Gyr old at 2 = 1, for 
the same cosmology as in Sect. 2. Hence it may hardly sus- 
tain a constant SFR larger than 10 2 Mq yr _1 for such long 
times. For our DSGs/DPSGs observed at z = 1.07, 1.39 or 
1.85, the maximum age of 1 Gyr corresponds to a maximum 
Zf of 1.39, 1.85 or 2.59, respectively. Hence these systems are 
mostly formed at 1 < 2 < 2. We consider additional model 
ages of 50, 100, 200, 300 and 500 Myr. 

PMOO consider also model starbursts with a SFR hav- 
ing an e-folding time of 100 Myr. These PMOO models are 
consistent with our models of starbursts with a short burst 
at an age of 100 Myr. However, PMOO do not consider an 
early post-starburst phase at 1 < 2 < 2. 

At this point the difference between a post-starburst 
galaxy and a passively evolving elliptical galaxy younger 
than 1 Gyr may sound a matter of semantics, since a DSG 
at 1 < 2 < 2 may evolve into an elliptical at 2 = and, even- 
tually, already at these redshifts. As far as stellar properties 
are concerned, the SEDs of the post-starburst stellar popu- 
lation under investigation and of an SSP with the same age 
and metallicity, identified as a passively evolving elliptical 
galaxy, are different (as expected) even after 1 Gyr from the 
onset of star formation, as shown in Fig. 1. However, it is the 
presence, distribution and structure of the (dusty) interstel- 
lar medium which provides fundamental physical differences 
between these two stellar systems in our models. 

The mass fractions of stars and residual gas are 74 (77) 
per cent and 22 (19) per cent, respectively, for a 1 Gyr- 
old stellar system formed in a short (long) burst. For DSGs 
with a short burst, these fractions are 87 per cent and 11 
per cent, respectively, at the end of the burst (i.e., at an 
age of 100 Myr). We assume that gas and dust components 
survive at least partly in model DSGs/DPSGs for the largest 
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Figure 1. The intrinsic SED of an SSP with Zq and age equal 
to 0.2 (red line) or 1 Gyr (magenta line) and the intrinsic SED of 
a stellar population originated in a 100 Myr burst with Zq and 
age equal to 0.2 (blue line), 0.3 (cyan line), 0.5 (green line) or 1 
Gyr (black line). Each SED is in the object's rest-frame and is 
normalized to the luminosity at 2.17 (im. It is evident that these 
two sets of SEDs are different at any intermediate age. 

time-window (1 Gyr) considered for these objects (see Sect. 
2.2.3). The fate of these components afterwards may channel 
the evolution of these starbursts towards either a "gas-poor" 
system (i.e., an early-type galaxy) or a "gas-rich" one (i.e., 
a late-type galaxy) at z — 0, both systems having a stellar 
mass of - 10 10 -10 n Mq. 

Thus a 1 Gyr-old DPSG at z = 1.85 could be the 
transition evolutionary phase of elliptical dusty formation 
occurred at intermediate redshift (z = 2-3) (Totani et 
al. 2001; Granato et al. 2001). However the dusty forma- 
tion of more massive Es at intermediate redshift would re- 
quire initial starbursts with SFR ^ 10 3 Mq yr _1 and/or 
hurst > 100 Myr (cf. Granato et al. 2001). 

2.2.2 Dust attenuation 

Gordon, Calzetti & Witt (1997 - hereafter referred to as 
GCW97) have investigated the nature of the dust in star- 
bursts by constructing different models of the stars and dust 
in DSGs, combining stellar evolutionary synthesis and ra- 
diative transfer through dust, and applying them to the ob- 
served UV/optical colours of 30 nearby starburst galaxies of 
the Calzetti sample (Calzetti, Kinney & Storchi-Bergmann 
1994). Their main conclusion are as follows. 

• Starburst dust has an extinction curve lacking the 
2175 A bump, like the Small Magellanic Cloud (SMC) curve 
and at variance with the Milky Way (MW) curve, and a 
steep far-UV rise, intermediate between these two curves. 
This is consistent with the conclusion of Calzetti et al.. 

• The dust/stars geometry that better explains the distri- 
bution of these 30 starbursts in various UV/optical colour- 



colour plots has an inner dust-free sphere of stars surrounded 
by an outer star-free shell with a two-phase clumpy, dusty 
medium. A geometry where dust and stars are uniformly 
mixed through the entire sphere can not explain the spread 
in the UV/optical colours, be its local distribution of dust 
homogeneous or clumpy. 

Models of dusty starbursts based on similar considerations 
well reproduce the optical colour-redshift distribution of 
Lyman-Break Galaxies at 2 < z < 4 (Vijh, Gordon & 
Witt 2003). In addition, Cimatti et al. (1997) suggest that 
the stellar component of ERO J164502+4626.4 should be 
strongly embedded in the dust in order to reproduce its ob- 
served extremely red colour. 

For all these reasons, we apply the Monte Carlo sim- 
ulations of the radiative transfer by WG00 for the SHELL 
geometry to DSGs/DPSGs at 1 < z < 2. This geometry 
assumes a dust/stars distribution where the stars extend to 
0.3 of the system radius and the dust extends from 0.3 to 1 
of the system radius. The dust in the shell is distributed in 
a two-phase clumpy medium, comprised of molecular clouds 
and diffuse gas, where the filling factor of the clumps is set 
to 0.15 and the low-density to high-density ratio is set to 
0.01. We consider the additional extreme hypothesis that 
dust is locally distributed in a single, homogeneous medium. 

We assume an SMC-like extinction curve and consider 
values of t v equal to 0.25, 0.5, 1, 2, 3, 4, 5, 6, 8, 10, 20, 30, 
40 and 50 for the two-phase clumpy dust distribution and 
equal to 0.25 through 8 for the homogeneous one. Here we 
recall that tv gives the total amount of dust in the shell, 
since the dust column density is proportional to tv x a x p, 
where a is the characteristic grain size and p is the density of 
the grain material. Values of tv at the low- and high-end of 
the distribution seem to be characteristic of starbursts and 
ULIRGs, respectively, at low- as well as high-z (e.g. Vijh et 
al. 2003 and references therein). 

We also assume that the gas emission at a given wave- 
length is attenuated by the same amount as the stellar emis- 
sion at that wavelength, whether the gas emission is either 
in a line or in the continuum. This is at odds with the re- 
sult that the stellar continuum suffers roughly half of the 
reddening suffered by the ionized gas (Calzetti et al. 1994; 
Fanelli, O'Connell & Thuan 1988; Mas-Hesse, Arnault & 
Kunth 1989) but the understanding of the attenuation of 
the ionized-gas line emission from dust is still controversial. 
However, the following results hold when no nebular emis- 
sion is taken into account, as demonstrated in App. A. 

Fig. 2 reproduces the attenuation curves (normalized to 
V band) for the SHELL dust/stars configuration and SMC- 
type dust locally distributed in a homogeneous (solid lines) 
or two-phase clumpy medium (dotted lines) as a function of 
tv- It also shows the values of the normalized attenuation 
function of nearby starbursts (Calzetti et al. 2000), known as 
"Calzetti law" , for the WG00 wavelengths (filled triangles) . 

Two important results are contained in Fig. 2. First, 
the attenuation function for the SHELL geometry with a 
two-phase clumpy distribution of SMC-type dust and tv = 
1.5 reproduces rather well the calibrated Calzetti law for 
nearby starbursts, as already shown by Vijh et al. (2003). 
Nevertheless a wide range of different attenuation curves 
may apply to individual DSGs/DPSGs, as a function of local 
distribution of the dust and dust content. We discuss the 
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effect of using the empirically determined Calzetti law for 
whatever dusty ERG (e.g. Cimatti et al. 2002a) in App. B. 

Second, the shape of the attenuation function depends 
on structure of the dusty medium and, in the two-phase 
clumpy case, on Ty. Changes in normalization also exist, 
in both cases, of course (not shown). In particular, a homo- 
geneous dusty medium provides a larger attenuation at any 
wavelength than a two-phase clumpy one. 

In addition, we stress that the relation between E(B — 
V)star and dust amount (i.e., ryj is usually quite non lin- 
ear and highly dependent on details of the structure of the 
dusty medium, as shown in Fig. 3. The existence of an initial 
"plateau" in E(B — V) s tar vs. ry in the two-phase clumpy 
case is due to fact that individual clumps become optically 
thick while the low-density interclump medium initially con- 
tributes little to the reddening. 

Finally, we note that, apart from the observationally- 
established attenuation function of Calzetti (1997) (see App. 
B), all the model attenuation functions considered by PM00 
do not contemplate dust with SMC-type extinction curve 
and distributed in a two-phase clumpy medium at the same 
time, as well as a SHELL-like geometry. We anticipate that 
these important differences do not affect our re-analisys of 
the two-colour classification criterion for ERGs. In fact our 
results are consistent with those of PM00 when continuously 
star-forming DSGs are considered, whatever structure of the 
dusty medium is assumed. 

2.2.3 The physics behind the SHELL configuration 

The existence of dust in young galaxies is expected because 
Type II supernovae (SNe II) are shown to produce dust 
grains (e.g., Dwek et al. 1983; Moseley et al. 1989; Kozasa, 
Hasegawa & Nomoto 1991; Todini & Ferrara 2001). Dust 
production from SNe II follows the formation and evolu- 
tion of massive stars on a nuclear timescale (few xlO 6 yr). 
Mass loss from evolved intermediate-mass stars contributes 
to dust formation as well (Gehrz 1989) and follows the 
slower evolutionary timescales of these stars (from ~ 10 8 
to ~ 10 9 yr). Dust is also destroyed by SN shocks (McKee 
1989; Jones, Tielens & Hollenbach 1996). Nevertheless, the 
detailed modelling of dust evolution in high-2: and nearby 
young starbursts shows that it is possible to build up a dust 
mass within a factor of two from total within a few x 10 7 yr 
(Hirashita & Ferrara 2002; Hirashita, Hunt & Ferrara 2002). 
Early dust is similar to the SMC one (Takeuchi et al. 2003) . 

Dust and residual gas face displacement and, eventually, 
partial/total removal from the system owing to the onset of 
galactic winds powered by SN explosions and accelerated 
up to velocities of several x 10 km s~ (e.g. Suchkov et 
al. 1994; Suchkov et al. 1996; Shu, Mo & Mao 2003). In 
the case of dust grains, destruction by sputtering must be 
taken into account (e.g. Ellison, Drury & Meyer 1997) . SN la 
will contribute to dust destruction/displacement as well as 
SN II, because their timescales span a wide range: roughly 
60 per cent of Type la SN progenitors have lifetimes ~ 3- 
4 x 10 s yr, while 30 per cent of them are expected to be 
very old (10 9 -10 10 yr) (Branch et al. 1995). The existence 
of galactic superwinds associated with dust extends from 
nearby to high-2: starbursts (Heckman 2001 and references 
therein; Adelberger et al. 2003). 

The fate of the surviving ejected dust grains is uncer- 




Figure 2. The attenuation function, normalized to its V-band 
value, (i.e., A x norm ) for: (a) the WTC92 CLOUDY geometry of 
a dusty elliptical galaxy with MW-type dust and Ty = 0.5 (long- 
dashed line); (b) the WG00 SHELL geometry of a DSG with 
SMC-type dust (solid or dotted lines for a homogeneous or two- 
phase clumpy dust distribution, respectively) and different values 
of t y . As a reference, we reproduce the normalized attenuation 
function for the WG00 SHELL model with clumpy SMC-type 
dust and Ty = 1.5 (short-dashed line), which is consistent with 
the normalized Calzetti law (filled triangles). 




two-phase clumpy _ 
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Figure 3. Stellar reddening E(B — V) s tar vs. ry for the SHELL 
geometry of a DSG with SMC-type dust distributed in a homo- 
geneous (solid line) or two-phase clumpy (dotted line) medium. 
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Figure 4. Distribution of model Es at z = 1.07 (a), z = 1.39 (b) 
or z = 1.85 (c) in the I c — K vs. K colour-magnitude diagram 
for the dust-free (thick solid line) or dusty (thin solid line) case. 
In each panel, the model age increases from left to right along 
each line. Long- and short-dashed lines reproduce the limiting 
magnitude of current ERO samples and the I c — K colour selection 
criterion for EROs, respectively. Note the extremely narrow locus 
populated by Es at 1 < z < 2 in this plane. 

tain for the following reasons. The ambient halo gas that 
the superwind shocks and the dense disk-gas that it en- 
trains travel at velocities of a few xlO 2 km s _1 , which are 
near to typical escape velocities (Nakai et al. 1987; Phillips 
1993). Furthermore, displaced grains may respond as much 
to radiation pressure (Davies et al. 1998; Norman & Ferrara 
1996) as they do to gas pressure (as well as to any poloidal 
magnetic field that might be present - Beck et al. 1994). 
Hence displaced dust may decouple from the gas so that its 
path through the halo becomes even more uncertain. Dis- 
placed grains may partially undergo both reconnection to 
the original disk and removal from the system, as well as 
they may exist in equilibrium positions above the plane of 
the galaxy (Davies et al. 1998; see also Greenberg et al. 
1987 and Barsella et al. 1989). This supports our assump- 
tion that the SHELL geometry applies also to the immediate 
post-burst phase of a DSG. 



3 PREREQUISITES FOR A MODEL BEING 
AN ERO 

Existing ground-based near-IR surveys (Drory et al. 2001; 
Cimatti et al. 2002b; Smail et al. 2002) reach a typical limit 
of K = 20.5. K < 20-21 is also the limit of validity of 
the PM00 method, introduced by these authors in order to 
avoid statistical contamination by z > 2 ellipticals in the 
optical/near-IR colour-colour region populated by DSGs. 

An important prerequisite of our models is to match 
observed photometric properties of EROs, i.e. K < 20.5 and 
e.g. 4 < I c — K < 8. Thus we show the distribution in the 



two-phase clumpy ISM homogeneous ISM 



6 
4 



B 

4 



6 
4 

2 



1 1 1 1 1 1 1 
- a 


i | i i i 

z-1.07 : 


b 


z=i.o7 : 


1 1 1 1 1 

c 

7 ' 


m m 

z-1.39 : 




: d .<^^ 

: ' ' ' I ' ' ' L 


I; I I 

z-1.39 : 




G 

r , , , 1 , , , 


Mill, 
z.1.85 : 


- f ,-j>? 


,1 .11 

/ \ 

z.1.85 : 



18 20 22 18 20 22 

K (mag) 



Figure 5. Distribution of model DSGs with a short burst and 
DPSGs at z = 1.07 (a, b), z = 1.39 (c, d) or z = 1.85 (e, f) in 
the I c — K vs. K colour-magnitude plane for the two-phase clumpy 
(a, c, e) and homogeneous (b, d, f) dust distributions. Long- and 
short-dashed lines are the same as in Fig. 4. Solid or dotted lines 
connect models with constant ry or age, respectively. In the upper 
left (right) corners of panels a,c,e (b,d,f) we reproduce the typical 
trend of models with increasing age (solid line) or increasing ry 
(dotted line) - see text. 
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Figure 6. The same as in Fig. 5 but for model DSGs with a long 
burst. Solid or dotted lines connect models with constant ry or 
age, respectively. In the upper left (right) corners of panels a,c,e 
(b,d,f) we reproduce the typical trend of models with increasing 
age (solid line) or increasing ry (dotted line) - see text. 
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I c — K vs. K colour-magnitude plane of our models of Es, 
DSGs with a short burst and DPSGs, and DSGs with a long 
burst in Figs. 4, 5 and 6, respectively. 

Here we remind that our model Es have a minimum stel- 
lar mass equal to 7.5 x 10 10 Mq whatever their age (from 
1 to 3-5 Gyr) is. Model DSGs with a short burst have stel- 
lar masses equal to 4.5 and 8.7 x 10 10 Mq when they are 
0.05 and 0.1 Gyr old, respectively, while model DPSGs have 
stellar masses equal to 8.1, 7.9, 7.7 and 7.4 x 10 10 Mq when 
they are 0.2, 0.3, 0.5 and 1 Gyr old, respectively. Finally, 
DSGs with a long burst have stellar masses equal to 0.4, 
0.9, 1.7, 2.5, 4.0 and 7.7 x 10 10 Mq when they are 0.05, 0.1, 
0.2, 0.3, 0.5 and 1 Gyr old, respectively. 

From Fig. 4 it is evident that model dust-free Es of 
7.5 x 10 10 Mq correspond to observed EROs only if they are 
at z < 1.4 (for any age at z ~ 1 or for ages less than 2 Gyr 
at z ~ 1.4). Dusty Es of the same stellar mass correspond 
to observed EROs only if they are at z ~ 1 (for any age). 
Hence stellar mass plays a key role: a ten times more massive 
passively evolving, old elliptical galaxy at 1 < z < 2 is 
always a candidate ERO, whether dusty or not. Indeed large 
near-IR surveys aimed at testing galaxy formation models 
are designed to detect Es more massive than ~ 10 11 Mq at 
high- 2 (e.g. Drory et al. 2001). 

The prerequisites for a model DSG/DPSG at 1 < z < 2 
being a candidate observed ERO are more complex, as illus- 
trated in Figs. 5 and 6 and summarised in Table 1. This is 
largely due to the well-known existence of the age/opacity 
degeneracy (Takagi, Arimoto & Vansevicius 1999 and refer- 
ences therein) . 

A dusty stellar system originated in a 100 Myr burst 
evolves towards redder I c — K colours when it ages, for fixed 
dust distribution and amount of dust. It also evolves towards 
brighter K magnitudes for the first 100 Myr and towards 
fainter K magnitudes during the following post-starburst 
phase. Conversely, it becomes fainter in K and redder in 
I c — K when ry increases, for a fixed age. In general, younger 
starbursts with tt urs t = 100 Myr must be more attenuated 
in order to match the properties of observed EROs, be their 
dusty medium two-phase clumpy or homogeneous (Fig. 5). 
In the second case, a lower amount of dust is needed in order 
to produce a given reddening (cf. Fig. 3) . 

A DSG with a 1 Gyr burst evolves towards brighter 
K magnitudes and redder I c — K colours when it ages, for 
fixed dust distribution and amount of dust. Conversely, it 
becomes fainter in K and redder in I c — K when tv increases, 
for a fixed age. Interestingly, starbursts with t burst = 1 Gyr 
do not match observed ERGs if their dusty medium is two- 
phase clumpy (in the way we assumed in Sect. 2.2.2) if not 
at z ~ 1 and for an age of 1 Gyr and ry ~ 30 (Fig. 6). 
They also hardly satisfy the ERO selection criterion at any 
age, if not for large amounts of dust (i.e., for ry > 20). Con- 
versely, DSGs with long bursts and a homogeneous dusty 
medium fall in the region populated by observed ERGs for 
low/intermediate amounts of dust and almost any age, at 
any redshift between 1 and 2. For these models the range 
in dustiness or age allowed by observations increases with 
decreasing z. We conclude that the model DSGs/DPSGs at 
1 < z < 2 span very well the region of the I c — K vs. K 
colour-magnitude plane populated by observed ERGs. 

At this point it is interesting to note what happens if a 
10 times lower mass (i.e., equal to 10 10 Mq) of primordial 



Table 1. Parameters of model DSGs/DPSGs selected as EROs 
with K < 20.5 



tburst 


age 


two-phase clumpy 


homogeneous 


z 


Myr 


Myr 








100 


50 


tv > 40 


3 < tv < 8 


1.07 


100 


100 


TY J> 40 


3 < tv < 8 


1.07 


100 


200 


20 < r v < 50 


2 < r v < 8 


1.07 


100 


300 


20 < tv < 40 


2 < rv < 8 


1.07 


100 


500 


10 < t v < 30 


2 < T V < fi 

£j < y V 


1.07 


100 


1000 


1 < tv < 20 


1 < TV < 6 


1.07 


100 


50 


t v > 30 


2 < r v ^ 6 


1.39 


100 


100 


30 < tv ^ 50 


2 < tv ^ 6 


1.39 


100 


200 




1 < tv < 6 


1.39 


100 


300 


2 < tv < 20 


1 < TV 5 


1.39 


100 


500 


1 < tv ^ 10 


1 ^ TV < 5 


1.39 


100 


1000 


0.25 ^ tv < 8 


0.25 <C tv ^ 4 


1.39 


100 


50 




2 < tv < 6 


1.85 


100 


100 




2 < tv < 6 


1.85 


100 


200 




1 < tv ^ 4 


1.85 


100 


300 


1 < TV ^ 5 


1 ^ tv < 4 
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0.25 < t v ^ 3 
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TV < 2 


tv < 2 
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3 $C i~\ t ?C 6 
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2 < tv < 8 


1.07 
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2 < tv < 8 


1.07 
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2 < tv < 8 


1.07 


1000 
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ry ~ 30 


2 < t v < 8 


1.07 


1000 


50 




2 < t v < 4 


1.39 


1000 
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2 < t v < 5 


1.39 


1000 


200 




2 < t v < 5 


1.39 


1000 


300 




2 < tv < 6 


1.39 


1000 


500 




2 < t v < 6 


1.39 


1000 


1000 




2 < tv < 6 


1.39 


1000 


300 




2 < tv < 3 


1.85 


1000 


500 




2 s£ tv < 3 


1.85 


1000 


1000 




1 < tv < 4 


1.85 



gas is transformed into stars at a constant rate on a timescale 
of 0.1 or 1 Gyr. In this case, we find that DSGs with a short 
(long) burst and SFR = 10 2 Mq yr^ 1 (10 Mq yr" 1 ) and 
their dusty post-starburst phases result to be too faint to be 
detected by present ground-based K-band surveys, whatever 
their redshift, dust amount and distribution are. Hence we 
conclude that these surveys are not sensitive to dusty for- 
mation of low-mass galaxies at 1 < z < 2 (if any). 



4 THE Ic K VS. J K COLOUR-COLOUR 
DIAGRAM 

Hereafter we will discuss only the I c — K vs. J — K colour- 
colour diagram since similar conclusions hold for the R c — K 
vs. J — K colour-colour diagram. 

The PM00 classification criterion for ERGs is: 

J-K = 0.36 (I c -K) + 0.46. 

ERGs with given observed I c — K colour are classified as Es 
(DSGs) at 1 < 2 < 2 if they have a bluer (redder) observed 
J — K colour than this. 
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Figure 7. The distribution in the I c — K vs. J — K colour-colour 
plane of dust-free (filled symbols) and dusty (empty symbols) 
Es at 1 < z < 2. In each panel model Es of selected ages are 
reproduced as triangles, squares and circles according to their 
rcdshifts. The short-dashed and solid lines represent, respectively, 
the ERO selection function and the PMOO classification criterion 
for ERGs. All model Es obey the latter criterion. 



4.1 Dust effects on the colours of passively 
evolving, old, dusty ellipticals 

Fig. 7 shows the distribution in the I c — K vs. J — K colour- 
colour plane of dust-free (filled symbols) and dusty (empty 
symbols) Es at 1 < z < 2. The short-dashed and solid lines 
represent, respectively, the ERO selection function (i.e., I c — 
K ^ 4) and the PMOO classification criterion for ERGs. 

All the model Es obey the PMOO classification criterion 
for ERGs, i.e., they fall on the left side of the delimiting 
line (solid line) in Fig. 7. This is re-assuring for dust-free 
Es, since the evolutionary models examined by PMOO adopt 
the same stellar tracks, while those adopted here (Maraston 
1998) rely on a different set of stellar models. 

The new result is that models of dusty Es still fall on 
the left of the PMOO separation line, for reasonable upper 
limits to the amount of dust present in these stellar systems. 

In general, the inclusion of dust attenuation moves the 
dust-free models in Fig. 7 towards redder I c — K and J — K 
colours, almost parallel to the age-sequence track, for fixed 
metallicity, at a given redshift, and closer to the PMOO sepa- 
ration line, whatever redshift and age of the model elliptical 
galaxy are. Thus dust attenution in Es may explain very red 
I c — K and J — K colours without the need of invoking ex- 
tremely old ages or higher redshifts (or higher metallicities 
- cf. PMOO). 

4.2 Dusty starburst/post-starburst galaxies 

Hereafter we refer to model DSGs with a 100 Myr burst and 
model DPSGs, both with a two-phase clumpy distribution 



of the dust, as basic model DSGs and basic model DPSGs, 
respectively, given the results of GCW97. 

4.2.1 DSGs with a short burst and DPSGs 

Fig. 8 reproduces the distribution of basic model 
DSGs/DPSGs with different age and dust amount in the 
I c — K vs. J — K colour-colour plane. It shows that DSGs 
with a 100 Myr burst at 1 < z < 2 obey the PMOO classifica- 
tion criterion, whatever the structure of their dusty medium 
and dust amount are. 

Conversely, an intermediate age, dusty post-starburst 
galaxy may be mistaken for either a passively evolving, old 
elliptical or a young, dusty starburst, both at 1 < z < 2, as 
a function of its redshift, when classified only on the basis 
of optical/near-IR colours such as I c — K (or R c — K) and 
J-K. 

In particular, basic model DPSGs with ages between 0.2 
and 1 Gyr populate the region (upper left quadrant in Fig. 

8) where only dust-free Es were expected to fall, according 
to PMOO (see also Mannucci et al. 2002). These "intruder" 
basic model DPSGs are at 1.3 < z < 2 and span a large 
range in tv (0.25 ^ ry ^ 40), according to their age. These 
values of ry are not unplausibly large, especially when the 
model DPSGs have a homogeneous dust distribution (Fig. 

9) . In this case, the reddening effect per given amount of 
dust is maximized (Fig. 3), so that these "intruders" have 
ry $s 5 (according to their redshift). We stress that the 
"intruder" DPSGs may be luminous enough to be observed 
in present surveys at K < 20.5, be their dust distribution 
homogeneous or two-phase clumpy (cf. Tab. 1). Hence the 
presence of a dusty post-starburst phase at 1.3 < z < 2 is a 
source of concern for the neat result of PMOO. 

Finally, we note that the mix of DPSGs with ages be- 
tween 0.2 and 1 Gyr at z ~ 1 and DSGs (with ages between 
0.05 and 0.1 Gyr) at 1 < z < 2 in Fig. 8 is due to a con- 
spiracy between redshift, shape of the intrinsic SED and be- 
havior of the attenuation function of these stellar systems. 
The class of DPSGs is a suitable candidate for interpreting 
the low dust emissions attributed to the bulk of dusty, bright 
EROs on the basis of submm observations (cf. Introduction) , 
as detailed in Sect. 5. 

4-2.2 DSGs with a long burst 

Figs. 10 and 11 show that model DSGs forming stars contin- 
uously in a long burst satisfy the PMOO classification crite- 
rion, be their dust distribution two-phase clumpy (but these 
models are too faint in K to be selected as EROs - see 
Fig. 6) or homogeneous, respectively. PMOO have considered 
analogous starbursts observed during their continuous star- 
formation activity. 

Hence we confirm that there is a neat separation be- 
tween passively evolving, old stellar populations (in a dust- 
free environment) and stellar populations formed in a con- 
tinuous way in a dusty starburst, both at 1 < z < 2, (PMOO) 
which holds whatever structure of the dusty medium and 
dust amount of DSGs are. 

We do not explicitly consider the case of early, dusty 
post-starburst phases of a DSG with a long burst since the 
results obtained for basic model DPSGs in Sect. 4.2.1 qual- 
itatively apply to this case. In particular, one may expect 
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Figure 8. The distribution of the basic model DSGs/DPSGs in 
the I c — K vs. J — K colour-colour plane, where lines are the same 
as in Fig. 7. In each panel models of selected age are reproduced as 
triangles, squares or circles according to their redshifts, with Ty 
increasing from 0.25 to 50 from the lower left to the upper right. 
Model DPSGs older than 200 Myr and with a large range in Ty 
not only populate the region of DSGs (upper right quadrant) but 
they also "intrude" the region (upper left quadrant) where only 
passively evolving, old, dust-free Es were expected to fall (PM00). 
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Figure 9. The distribution of DSGs with a short burst /DPSGs 
with a homogeneous dusty medium in the same plane as in Fig. 
8. Each panel shows models of selected age at different z, with 
tv increasing from 0.25 to 8 from the lower left to the upper 
right. The "intruder" (see text) DPSGs with a homogeneous dusty 
medium have stellar populations similar to those of basic model 
DPSGs in Fig. 8. 



Figure 10. The distribution of model DSGs with a 1 Gyr burst 
and a two-phase clumpy dust distribution in the I c — K vs. J — K 
colour-colour plane as a function of age and z. In each panel Ty 
increases from 0.25 to 50 from the lower left to the upper right. 
Lines and symbols are the same as in Fig. 8. These DSGs meet 
the PM00 classification criterion for ERGs, even though they are 
too faint to be selected as EROs (Fig. 6). 
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Figure 11. The same as in Fig. 10 but for a homogeneous dust 
distribution. Each panel shows models of selected age at different 
z, with ry increasing from 0.25 to 8 from the lower left to the 
upper right. DSGs with a 1 Gyr burst and a homogeneous dusty 
medium meet the PM00 classification criterion for ERGs. These 
DSGs may be bright enough to be selected as EROs (cf. Tab. 1). 
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that this new class of DPSGs presents "intruders" which 
have a lower dust amount than "intruder" basic model 
DPSGs observed at the same time after the end of the 
burst and at the same redshift, owing to the existence of 
the age/opacity degeneracy. This is indeed what we find. 
As an example, a 1.2 Gyr-old dusty post-starburst phase 
following a burst with constant SFR = 10 2 Mq yr" 1 on a 
timescale of 1 Gyr and observed at z — 1.39 has I c — K = 4.44 
(5.10) and J — K = 1.95 (2.31), if its clumpy dusty medium 
has tv — 2 (20). This object falls well into the E domain 
(PM00) for tv = 2 and, in general, is an "intruder" for 
low/intermediate values of tv; it falls into the DSG domain 
(PM00) for high values of tv (e.g. 20 or less so). Conversely, 
a 0.3 Gyr-old basic model DPSG observed at z = 1.39 is an 
"intruder" when its clumpy dusty medium has 2 < tv < 20 
(cf. Tab. 1). 

4.3 The case of sub-solar metallicities 

From our further modelling and from the analysis of PM00 
(their Fig. 3c), it is easy to realise that the previous conclu- 
sions extend to Es and DSGs/DPSGs with metallicities dif- 
ferent from solar. For example, we find "intruders" among 
DPSGs with a two-phase clumpy distribution of the dust 
and a tenth-solar metallicity (not shown). In this case, the 
"intruders" are older than 300 Myr and have values of tv 
between 0.50 and 30 (according to their z). 

4.4 The SED of an "intruder" DPSG 

Here we give an example of the non-negligible subsample 
of computed models of Es and DPSGs which populate the 
same region in the I c — K vs. J — K colour-colour plane (cf. 
Figs. 7, 8 and 9). Fig. 12 shows the normalized SED of an 
SSP (identified as a passively evolving elliptical galaxy) with 
Zq and age equal to 4 Gyr at z = 1.39 (thick solid line) and 
the normalized SED of a basic model DPSG with Zq, age 
equal to 500 Myr and tv = 6 at z = 1.85 (thin solid line). 
Each SED is in the observed frame and is normalized to the 
luminosity of the object at 2.17 /jm. The spectral coverages 
of the I c , J and K filters are shown as dashed areas. The 
visual inspection of Fig. 12 shows that these two objects 
have identical I c — K and J — K colours. Indeed I c — K = 5.01 
and J — K = 2.11 for the elliptical galaxy and I c — K = 4.97 
and J - K = 2.09 for the DPSG 2 . 

In theory, complementary optical photometry in U, B 
and V may help distinguishing DPSGs at 1.3 < z < 2 from 
Es at 1 < z < 2, which can not be distinguished by the PM00 
classification criterion. Model Es have e.g. values of V — R c 
ranging from 2 to 3 when their V — K ranges from 7.5 to 
12, the former range being roughly 1 mag redder than that 
of model DPSGs with the same V — K colour. However, our 
models of the oldest Es and of the oldest, dustiest DPSGs, 
both at 1 < 2 < 2, predict apparent magnitudes fainter than 
28 V-mag, so that only the brightest (21 < V < 26) among 

2 The elliptical galaxy has K = 20.73 (for a stellar mass of 

7.5 X 10 10 Mq), while the DPSG has K = 21.04. Hence these 
two models correspond to objects beyond the limiting K-band 
magnitude of the present ERO samples. N.B.: this is only a mat- 
ter of normalization. 
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Figure 12. The SED of an SSP (i.c, an elliptical galaxy) with 
Zq and age of 4 Gyr at z = 1.39 (thick solid line) and the SED 
of a basic model DPSG with Zq, age of 500 Myr and Ty = 6 at 
z = 1.85 (thin solid line). Each SED is in the observed frame and 
is normalized to the luminosity at 2.17 (im. The spectral coverages 
of the I c , J and K filters are shown as dashed areas. It is evident 
that these two objects have similar I c — K and J — K colours. 



these objects (either non-extremely attenuated DPSGs at 
1 < z < 2 or very massive Es at z ~ 1), seem to be easily 
detectable, given the current observational capabilities. 

Hence we conclude and confirm that well-sampled SEDs 
are needed for a robust classification of ERGs (e.g. Smith et 
al. 2001; Smail et al. 2002; Miyazaki et al. 2002), as well as 
observations in X-rays, radio and/or far-IR/submm (see In- 
troduction). Photometry in the far-IR/submm is suggested 
almost by default, even though present observational cam- 
paigns with SCUBA (Andreani et al. 1999; Mohan et al. 
2002) show that the bulk of bright (K < 19.5-20) ERGs has 
total mid-IR-mm rest-frame luminosities lower than those 
of ULIRGs. We discuss this in Sect. 5. 

Finally, a way of distinguishing DPSGs from Es, both 
at high-z, may be provided by the presence/absence of 
intermediate-age (i.e., 10 8 -10 9 yr-old) stellar populations. 
We will discuss this in a future paper. 



4.5 On the classification of ERGs 

The existence of dusty ERGs mixed with dust-free objects 
having passively evolving, old stellar populations (i.e., Es) 
and vice versa has been recently claimed by Miyazaki et 
al. (2002), following the application of the photometric red- 
shift technique to multi- wavelength photometry obtained for 
247 EROs (see Fig. 5 in their paper). Evidence in favour of 
the diversity of ERGs comes also from the multi- wavelength 
analysis of 68 EROs by Smail et al. (2002) via a similar 
photometric redshift technique. Nevertheless, these authors 
conclude that the PM00 classification criterion is reasonably 
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effective, at least in the R c — K vs. J — K colour-colour plane 
(see Fig. 5 in Smail et al.). 

Fig. 13 reproduces the distribution in the I c — K vs. 
J — K colour-colour plane of 45 ERGs out of the ERO sam- 
ple of the previous authors. This subsample comprises ERGs 
brighter than K = 20.5 (5a) and detected at least in the I c 
or J band. ERGs classified by Smail et al. as evolved (i.e., 
with passively evolving, old stellar populations in a dust- 
free environment), dusty (i.e., starburst) or other (i.e., with 
SEDs not matching those of the previous two classes) are re- 
produced as empty circles, filled circles and empty triangles, 
respectively. The visual inspection of Fig. 13 shows that: 

• half of the dusty objects fall within la from the PM00 
separation line (on both sides). The existence of potential 
confusion near this separation line was stressed by PM00. 

• Two objects out of the 12 ERGs classified by Smail et 
al. as galaxies with evolved SEDs fall within la from the 
PM00 separation line (on both sides) . 

• Two more ERGs with evolved SEDs fall on the right 
side of the PM00 separation line. They could be evolved 
galaxies at z > 2 (see PM00). 

• ERGs with SEDs more complex than evolved- or dusty- 
type fall everywhere. 

For passively evolving, old ellipticals and dusty starbursts at 
1 < z < 2, the PM00 classification criterion, confirmed by 
our models (see Sect. 4.1, 4.2 and 4.3), is statistically effec- 
tive, as pointed out by Smail et al. (2002). Indeed also the 
median colours of these two galaxy populations are separate 
(see Table 2 in Smail et al.). 

The more complex distribution in Fig. 13 (see also 
Miyazaki et al. 2002) is only partly accounted for by our 
theoretical exploration. In fact we find that model DPSGs 
may overlap with model Es, independent of their dust dis- 
tribution (Sect. 4.2.1) and metallicity (Sect. 4.3). 

Peculiar systems (e.g. AGN-host galaxies, mergers) or 
bulge+disk galaxies (with different composite stellar pop- 
ulations, affected by dust in different ways) at 1 < z < 2 
may also be part of the ERO population. In fact the latter 
objects have been found by Smith et al. (2002), Miyazaki et 
al. (2002) and Yan & Thompson (2003). 

The set of models discussed here does not account for 
these dusty galaxies. However we note that a passively evolv- 
ing, old, dust-free population associated with the bulge com- 
ponent of a bulge+disk galaxy may present redder colours 
than the same population associated with a dust-free el- 
liptical galaxy, owing to the attenuation from the dust 
in the disk. At the same time, the disk stellar popula- 
tion may behave like a continuously star-forming DSG with 
low/intermediate SFRs (cf. PM00). Hence a bulge+disk 
galaxy at 1 < z < 2 may be expected to fall on both sides 
of the PM00 delimiting line. We defer to a future paper the 
investigation of normal bulge+disk galaxies at z > 1 as can- 
didate ERGs. 



4.6 On the photometric definition of EROs 

Fig. 13 shows that 3 objects (with different classifications), 
selected as EROs by Smail et al. (2002) on the basis of their 
R c — K colour, do not match the I c — K colour selection 
criterion of EROs at more than la, as already noted by 
these authors. Smail et al. find that their field contains 68 



(from Smail et al. 2002) 




J-K (mag) 

Figure 13. The distribution of 45 ERGs observed by Smail et al. 
(2002) in the I c — K vs. J — K colour-colour plane. ERGs classified 
by Smail et al. as evolved, dusty or other are reproduced as empty 
circles, filled circles and empty triangles, respectively. The short- 
dashed and solid lines represent the ERO selection function and 
the PM00 classification criterion for ERGs, respectively. Note that 
these objects were identified as EROs by Smail et al. on the basis 
of their R c — K colour. 



EROs with R c - K > 5.3 against 99 objects with I c - K > 4 
and suggest that the latter is a less stringent definition of 
an unusually red source. In fact, 63 (93 per cent) of these 68 
R c — K colour-selected EROs have I c — K > 4. 

With our models, we can investigate the nature of these 
"ambiguous" objects. We find that model Es at 1 < z < 
2, brighter than K = 20.5, obey both the R c - K > 5.3 
and I c — K ^ 4 colour selection criteria for EROs, whether 
dusty or not, for both solar and half-solar metallicities (not 
shown). Conversely, model DSGs/DPSGs at 1 < z < 2, 
brighter than K = 20.5, selected via either I c — K or R c — K, 
produce synthetic ERO samples with different objects and 
sizes, for both solar and a tenth-solar metallicities. This is 
shown in Figs. 14 and 15 for Zq DSGs with a short burst and 
DPSGs, both with a two-phase clumpy dust distribution, 
and for Zq DSGs with a long burst and a homogeneous 
dust distribution, respectively. 

Interestingly, the R c — K ^ 5.3 colour selection criterion 
systematically excludes more basic model DSGs/DPSGs at 
1.3 < 2 < 2 (mostly 300 Myr-old DPSGs with 0.50 < 
tv < 10, but also DSGs with -ry = 40). Conversely, the 
I c — K ^ 4 colour selection criterion tends to exclude 
more DSGs/DPSGs at z ~ 1 (mostly 500 Myr-old DPSGs, 
with 4 tv ^ 20). In case of model DSGs with a short 
burst and DPSGs with a homogeneous dust distribution, 
the R c — K ^ 5.3 colour selection criterion systematically 
excludes 300 Myr-old DPSGs with t v = 1 and 500 Myr-old 
DPSGs with tv — 0.50, both at z > 1.8. Conversely, there is 
no evidence for a bias in z, caused by the different colour se- 
lection criteria, for model DSGs with a long burst, whatever 
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the structure of their dusty medium is (not shown). How- 
ever, these DSGs are excluded by the R c — K ^ 5.3 colour 
selection criterion if they have ry = 30 (two-phase clumpy 
case) or 2 ^ ^ 3 (homogeneous case). Thus we conclude 
that selection via the R c — K colour bias the resultant ERO 
sample towards a larger ratio of Es (at any z between 1 and 
2) to DSGs/DPSGs at 1.3 < z < 2; vice versa, selection via 
the I c — K colour bias it towards a larger ratio of Es (at any 
z between 1 and 2) to DSGs/DPSGs at z ~ 1. 

In addition, we find that ERO selection via the R c — K 
colour tends to exclude DSGs/DPSGs brighter than K = 
20.5 at 1.3 < z < 2, if they have relatively low amounts of 
dust (i.e., with 0.50 ^ iy ^ 3), whatever the structure of 
their dusty medium is. 

Hence large discrepancies may be expected in the total 
number of selected EROs and in the statistical properties 
of the different sub-populations of ERGs present in a given 
ERO sample, whether selection is via R c — K or I c — K. These 
considerations help understanding the 33 per cent mismatch 
in the number of EROs selected via either R c — K or I c — K 
found by Smail et al. (2002). 

We note that Yan & Thompson (2003) conclude that 
the Es-to-DSGs ratio in a given ERO sample is larger, when 
selection is via the R c — K colour instead of the I c — K 
colour, on the basis of the observed optical/near-IR colours 
expected for model stellar populations (Bruzual & Chariot 3 ) 
approximating a passively evolving, old elliptical or simulat- 
ing galaxies of different ages with extended star formation 
activity, all at z ~ 1. They do not consider dust attenu- 
ation, however. Our results support and complement their 
conclusion. 



5 TOTAL MID-IR-(SUB)MM EMISSION OF 
DUSTY STARBURST/POST-STARBURST 
GALAXIES AT 1< Z < 2 

Here we detail on the relation between amount of dust and 
structure of the dusty medium, on one side, and total mid- 
IR-(sub)mm emission from dust, on the other side, for model 
DSGs/DPSGs. The radiative transfer models described in 
Sect. 2.2.2 do not model re-emission from dust, and, thus, 
we can not produce and analyse SEDs in the rest-frame 
mid-IR-(sub)mm for these stellar systems. However, they 
provide a robust estimate of the total emission of model 
DSGs/DPSGs at these wavelengths, as given by the amount 
of light absorbed by internal dust in the wavelength range 
0.1-3 /im (in the rest-frame). We refer to this estimate as 

Ldust ■ 

We find that Ldust has a dynamic range of more than 
500 (from ~ 10 10 to ~ 5.6 x 10 12 L ) for DPSGs through 
DSGs, as a function of the dust parameter space. In partic- 
ular, Ldust ranges between 2 and 5.6 x 10 12 Lq for DSGs 
with a short burst, between 2 and 7.6 x 10 1 Lq for DSGs 
with a long burst, and between 10 10 and 5.7 x 10 11 Lq for 
DPSGs. For fixed dust parameters, Ldust increases (up to 50 
per cent) with increasing age for both DSGs, but decreases 
(down to an order of magnitude) with increasing age for 
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Figure 14. The distribution in the R c — K vs. I c — K colour— colour 
plane of dusty stellar systems formed in a short burst, of Zq and 
different ages, at 1 < z < 2, brighter than K = 20.5 and with a 
two-phase clumpy dust distribution. The two short-dashed lines 
represent the R c — K ^ 5.3 and I c — K ^ 4 selection functions of 
EROs. ERO selection via either colour criterion produces samples 
with different objects and sizes. A selection effect in z is present. 



long burst, homogeneous ISM 
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Figure 15. The distribution in the R c — K vs. I c — K colour— colour 
plane of model DSG with a long burst, of Zq and different ages, 
at 1 < z < 2, brighter than K = 20.5 and with a homogeneous 
dust distribution. Conclusions are the same as in Fig. 14, but for 
the absence of a selection effect in z. 
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DPSGs. Conversely, for fixed age, Ldust increases by a fac- 
tor of about three when ry increases from 0.25 to 50 (from 
0.25 to 8) for the two-phase clumpy (homogeneous) case, for 
both DSGs and DPSGs. 

The reason for the high dynamic range in Ldust is 
twofold. On one hand, a dusty medium with given dust type 
absorbs stellar light at a given wavelength with an efficiency 
which strongly depends on dust amount, structure of the 
dusty medium, and distribution of this medium with respect 
to the stellar sources. Hence all these aspects contribute to 
determining the re-emission from dust grains in the mid-IR- 
(sub)mm wavelength domain. This is investigated in details 
by Misselt et al. (2001) for the same dust/stars geometry 
of a DSG/DPSG as that assumed in Sect. 2.2.2. Here it is 
important to stress that model DSGs/DPSGs selected as 
ERGs have very different dust properties (cf. Table l) 4 . 

On the other hand, dust absorbs and scatters more ef- 
ficiently UV photons than optical/near-IR ones. Available 
optical/near-IR photons may replace available UV photons 
as the major source of dust heating, if the optical depth at 
optical/near-IR wavelengths of the dusty medium is large 
enough. For a model DPSG, the intrinsic UV-to-optical lu- 
minosity ratio decreases towards when the model age ap- 
proaches 1 Gyr (Fig. 1). Thus ageing DPSGs will tend to 
have lower and lower total mid-IR-(sub)mm luminosities, 
for fixed dust parameters. 

For the dust parameter space explored here, only model 
DSGs with SFR = 10 3 M yr _1 over a timescale of 100 
Myr may be associated with Ldust ^ 10 12 Lq and, thus, 
be classified as ULIRGs. Conversely, model DSGs with 
SFR = 10 2 M yr 1 over a timescale of 1 Gyr typically 
have total mid-IR-(sub)mm luminosities of a few to sev- 
eral xlO 11 Lq. DPSGs may be associated with even lower 
values of dust emission (down to ~ 10 10 Lq), especially 
when their post-burst stellar population approaches 1 Gyr. 
Hence we conclude that the bulk of model DSGs/DPSGs 
does not belong to the ULIRG class. This is true also for 
model DSGs/DPSGs selected as EROs with K < 20.5 (cf. 
Table 1) or brighter (cf. Figs. 5 and 6). 

Our conclusion is consistent with those of Andreani et 
al. (1999) and Mohan et al. (2002), which are based on ob- 
servational results for a slightly brighter ERO sample (see 
Introduction) . 



6 DISCUSSION AND CONCLUSIONS 

In the previous sections we have modelled the spec- 
tral energy distributions (SEDs) of simple stellar popu- 
lations, intended to simulate passively evolving, old (> 
1 Gyr) elliptical galaxies (Es), and the SEDs of young 
(^ 0.1 Gyr) /intermediate-age (0.1-1 Gyr), dusty star- 
burst galaxies (DSGs) and of intermediate-age, dusty post- 
starburst galaxies (DPSGs), all at 1 < z < 2, as candidate 
extremely red objects (EROs). 

In the early post-starburst phases, massive stars are 
absent (owing to the absence of star formation) but 

4 In particular, we note that structure of the dusty medium may 
be traded for dust amount in producing a given value of Ldust- 
This is a caveat against the simplistic use of dust emission as a 
proxy of dust amount in dusty ERGs. 



intermediate-age stars are present. For these phases, we as- 
sume that at least part of the dusty medium of the previous 
starburst phase survives during a maximal time-window of 
0.9 Gyr. 

We have investigated observed optical/near-IR colours 
and magnitudes of the previous models and, in particular, 
their distribution in the I c — K vs. J — K colour-colour 
plane. In this colour-colour diagram, Pozzetti & Mannucci 
(2000) claim that passively evolving, old stellar popula- 
tions in a dust-free environment (i.e., elliptical galaxies) and 
stellar populations produced continuously in a dusty star- 
burst galaxy are reasonably well separated, whatever the 
dust/stars configuration is. 

This study complements the seminal analysis of 
Pozzetti & Mannucci for four main reasons. 

• First, we make the new case of dusty post-starburst 
galaxies at 1 < z < 2 as candidate EROs. In particular, wc 
discuss the early, dusty post-starburst phases of stellar sys- 
tems formed (mostly) at 1 < z < 2, in a single burst of finite 
duration (tburst) of the order of 0.1 Gyr (short burst). The 
results obtained for these DPSGs may be easily extended to 
coeval, early, dusty post-starburst phases of galaxy forma- 
tion in a single burst with e.g. hurst = 1 Gyr (long burst), 
when the existence of a degeneracy between age and opacity 
of the models is taken into account. 

We note that the DSGs with long or short bursts under 
investigation sustain star formation rates (SFRs) of the or- 
der of 10 2 M yr -1 or 10 3 M Q yr -1 , respectively. SFRs 
~ 10 3 Mq yr -1 have been inferred from observations of 
high-z (sub)mm sources (e.g., Smail, Ivison & Blain 1997; 
Hughes et al. 1998; Barger et al. 1998; Cimatti et al. 1998; 
Eales et al. 1999; Bertoldi et al. 2000; Lutz et al. 2001; Smail 
et al. 2003). 

The presence of a considerable population of starburst 
galaxies with large SFRs over a short (~ 100 Myr) timescale 
at 1 < z < 2 is qualitatively consistent with the in- 
terpretation of the Lick absorption-line indices of nearby 
intermediate-mass (i.e., with a velocity dispersion of ~ 
200 km s _1 ) early- type galaxies by Thomas et al. (2002). 

• Second, we explore the impact of dust attenuation 
(Witt et al. 1992) on the optical/near-IR colours of dusty 
Es at 1 < 2 < 2. 

• Third, for all DSGs/DPSGs, we adopt a new dust/stars 
geometry (Witt & Gordon 2000) and explore the impact of 
dust distribution (two-phase clumpy vs. homogeneous) and 
dust amount on the observed optical/near-IR colours. This 
geometry is consistent with observed properties of nearby 
starbursts and Lyman Break Galaxies at 2 < z < 4. 

• Fourth, we calculate the total emission from dust in the 
rest-frame mid-IR-(sub)mm for DSGs/DPSGs. 

The distribution of model Es and DSGs/DPSGs at 1 < 
z < 2 in the optical/near-IR colour-magnitude and colour- 
colour diagrams well reproduces that of observed extremely 
red galaxies (ERGs) (Sect. 3 and 4). Hereafter we summarise 
our results. 

• We find that DPSGs, with ages between 0.2 and 1 Gyr, 
at 1.3 < z < 2 may mix with Es at 1 < z < 2 in the 
I c - K vs. J — K colour-colour diagram, for a large range in 
dust amount. This holds whatever the structure (two-phase 
clumpy or homogeneous) of their dusty medium is. This re- 
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suit is a source of concern for the two-colour classification 
criterion of Pozzetti & Mannucci. 

• Conversely, intermediate-age DPSGs at z ~ 1 may mix 
with DSGs of age less than 1 Gyr at 1 < z < 2 in the same 
colour-colour plot, for a large range in dust amount and 
whatever the structure of the dusty medium is. This result 
and the previous one could explain in part the complexity 
of the ERG classification recently emerged from the analy- 
sis of statistical samples of EROs, where this classification 
was based on multi-wavelength photometry (e.g. Smail et 
al. 2002; Miyazaki et al. 2002) and not on optical/near-IR 
colour-colour plots. 

• On the other hand, we confirm the result of Pozzetti & 
Mannucci. In addition, we find that this result holds whether 
the dusty medium of a DSG is two-phase clumpy or homo- 
geneous and whether high-z ellipticals are dust-free or as 
dusty as nearby ones. 

• Model Es at 1 < z < 2 populate a narrow locus in the 
I c — K vs. J — K colour-colour plane, whether dusty ot not. 
This means that dust attenution may explain very red I c — K 
and J— K colours for ellipticals at 1 < z < 2 without the need 
of invoking extremely old ages and/or higher metallicities 
and/or higher redshifts. 

• The picture is the same when the distribution of the 
models in the R c — K vs. J — K colour-colour plane is inves- 
tigated (cf. Pozzetti & Mannucci 2000). Hence we conclude 
and confirm that well-sampled SEDs are needed for a robust 
classification of ERGs (e.g. Smith et al. 2001; Smail et al. 
2002; Miyazaki et al. 2002), as well as observations in X- 
rays (Alexander et al. 2002; Brusa et al. 2002; Smail et al. 
2002; but see Stevens et al. 2003), radio (Smail et al. 2002) 
and/or far-IR/submm (Cimatti et al. 1998; Andreani et al. 
1999; Mohan et al. 2002). 

• We find that the total rest-frame mid-IR-(sub)mm lu- 
minosity ranges from ~ 10 10 to ~ 5.5 x 10 12 L for DPSGs 
through DSGs. This high dynamic range is consistent with 
the large range in total rest-frame mid-IR-(sub)mm lumi- 
nosity estimated for ERGs observed in the submm (Cimatti 
et al. 1998; Andreani et al. 1999; Mohan et al. 2002). 

This study offers two viable, not mutually exclusive in- 
terpretations to the observational results and conclusions of 
the previous authors. One interpretation is that re-emission 
by dust grains in the mid-IR-(sub)mm wavelength domain 
depends on dust type and amount, structure of the dusty 
medium and dust/stars configuration, since the absorption 
efficiency of stellar light at given wavelength depends on 
these parameters. We find that the extremely red colours 
of DSGs/DPSGs at 1 < z < 2 do not necessarily indicate 
large amounts of dust (Sect. 3), as instead believed for e.g. 
SCUBA submm sources, and may be reproduced under a 
very large range in structure of the dusty medium. Hence 
the submm fluxes of DSGs/DPSGs at 1 < z < 2 are not ex- 
pected to be necessarily similar to those of SCUBA submm 
sources at similar/higher z. 

The other interpretation is the reduced contribution of 
UV photons to dust heating in the model DPSGs, which 
seem to constitute a non negligible, candidate ERO sub- 
population. Dust absorbs and scatters more efficiently UV 
photons than optical/near-IR ones. Available optical/near- 
IR photons may replace available UV photons as the major 
source of dust heating if the optical depth at optical/near- 
IR wavelengths of the dusty medium is large enough. For 



a model DPSG, the intrinsic UV-to-optical luminosity ratio 
decreases towards when its age approaches 1 Gyr (Fig. 1). 
Thus ageing DPSGs will tend to have lower and lower to- 
tal rest-frame mid-IR-(sub)mm luminosities, for fixed dust 
parameters. 

• Furthermore, for the dust-parameter space consid- 
ered here, we find that only model DSGs with SFR = 
10 3 Mq yr" 1 over a timescale of 100 Myr may be associated 
with a total mid-IR-(sub)mm luminosity (in the rest frame) 
greater than 10 12 Lq and, thus, be classified as ULIRGs. 
Conversely, model DSGs with SFR = 10 2 M yr _1 over 
a timescale of 1 Gyr typically have total mid-IR-(sub)mm 
luminosities of a few to several xlO 11 Lq. DPSGs may be 
associated with even lower values of dust emission (down 
to ~ 10 1() Lq). We conclude that the bulk of dusty, bright 
(K < 20.5) ERGs is not associated with ULIRGs, consistent 
with the conclusion of Andreani et al. (1999) and Mohan 
et al. (2002), based on observations of a slightly brighter 
sample. 

• Typical ground-based surveys with K < 20.5 (e.g., 
Drory et al. 2001; Cimatti et al. 2002b; Smail et al. 2002) 
are biased against the following objects at 1 < z < 2 (if 
any): ellipticals (dusty or not) with stellar masses lower 
than ~ 10 11 Mq; dusty starbursts with e.g. a constant 
SFR ^ 10 2 Mq yr" 1 (10 Mq yr" 1 ) on a timescale of 0.1 
(1) Gyr, identified as events of low-mass galaxy formation 
at these redshifts; early, dusty post-starburst phases of these 
dusty starbursts. 

• Model Es (dusty or not) at 1 < z < 2 are always se- 
lected as EROs (if massive enough) , whether selection is via 
the I c — K or R c — K colour. 

• Conversely, selection of observable (i.e., with K < 20.5), 
coeval DSGs/DPSGs as EROs is a complex function of the 
properties of their dusty medium. The explanation is that 
the extremely red colours (R c — K > 5.3 or I c — K ^ 4) 
of model DSGs/DPSGs at 1 < z < 2 derive from radiative 
transfer through their dusty media. 

• We find that selection via the R c — K colour bias the 
resultant sample of observable EROs towards a larger ra- 
tio of Es (at any z between 1 and 2) to DSGs/DPSGs at 
1.3 < z < 2; vice versa, selection via the I c — K colour bias 
it towards a larger ratio of Es (at any z between 1 and 2) 
to DSGs/DPSGs at z ~ 1. This result confirms and comple- 
ments that of Yan & Thompson (2003). 

Hence important discrepancies may be expected in the 
total number of selected EROs and in the statistical prop- 
erties of the different sub-populations of ERGs present in 
a given ERO sample, whether selection is via R c — K or 
Ic - K. For example, DSGs/DPSGs at 1.3 < z < 2 with 
the lowest amounts of dust are not selected as EROs via the 
R c — K ^ 5.3 colour criterion. 

These considerations help understanding the 33 per cent 
mismatch in the number of EROs with K < 20.5 selected 
via either R c — K or I c — K found by Smail et al. (2002) . 
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Figure Al. I c — K vs. J — K for Zq model DSGs with a short 
burst and DPSGs, all of them with a two-phase clumpy dusty 
medium, as a function of age and z, when no nebular emission is 
taken into account, as well as its attenuation from dust. In each 
panel ry increases from 0.25 to 50 from the lower left to the upper 
right. Lines and symbols are the same as in Fig. 8. Conclusions 
are the same as those reached from Fig. 8. 
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APPENDIX A: THE NON-IMPORTANCE OF 
THE NEBULAR CONTINUUM 

We reproduce the distribution in the I c — K vs. J — K colour- 
colour plane of Zq DSGs with a short burst and DPSGs 
(Fig. Al) and DSGs with a long burst (Fig. A2), all of them 
with a two-phase clumpy dust distribution, as a function 
of age and redshift, when no nebular emission is computed. 
Hence no attenuation of the gas emission, either in the line 
or in the continuum, from dust is computed as well. Mod- 
els are represented as in Figs. 8 and 10, respectively, where 
nebular emission was taken into account together with its 
attenuation from dust as described in Sect. 2.2.2. 

From the comparison of these four figures it emerges 
that nebular emission and its attenuation do not play a role 
in the selection of a DSG/DPSG at 1 < z < 2 as an ERO. 
They do not affect the presence of "intruder" DPSGs either. 
Hence we conclude that the main results in Sect. 4.2 do 
not depend on our choice for the attenuation of the nebular 
emission from dust. 



Figure A2. The same as in Fig. Al but for DSGs with a long 
burst. In each panel ry increases from 0.25 to 50 from the lower 
left to the upper right. Conclusions are the same as those reached 
from Fig. 10. 



APPENDIX B: ON THE APPLICATION OF 
THE CALZETTI LAW TO DUSTY ERGS 

Calzetti et al. (1994) have empirically constrained an extinc- 
tion function from the UV and optical spectra of 39 starburst 
and blue compact galaxies at z = 0. This extinction curve 
has an overall UV/optical slope flatter than the MW one 
and lacks the 2175 A dust absorption feature, as confirmed 
by GCW97. Further inclusion of IR broad-band data for 19 
of these 39 galaxies helped constrain the shape of the char- 
acteristic (i.e., average) attenuation function, but not fix its 
zero point (Calzetti 1997). The absolute calibration of the 
attenuation function of DSG at z — was possible owing to 
the availability of far-IR photometry for a subsample of only 
5 objects (Calzetti et al. 2000). In its final layout, the so- 
called "Calzetti law" has a fixed wavelength dependence (as 
shown in Fig. 2) and is normalized to the stellar-continuum 
reddening E(B — V) s tar- Since its discovery, the Calzetti law 
has been used to describe dust attenuation in DSGs at high- 
z (e.g., Cimatti et al. 2002a) under the intrinsic assumption 
that dust properties and content, as well as dust/stars dis- 
tribution and structure of the dusty medium, are similar in 
high- z DSGs and in the nearby starbursts of the Calzetti 
sample. 

Thus we produce new models of Zq DSGs with a short 
burst and DPSGs (Fig. Bl), and DSGs with a long burst 
(Fig. B2) by adopting the Calzetti law instead of the simu- 
lated attenuation curves described in Sect. 2.2.2. Two vales 
of the stellar-continuum reddening (i.e., 0.5 and 0.8) are 
considered, as in PM00 and consistent with Cimatti et al. 
(2002a). Even in this case we find that model DPSGs par- 
tially "intrude" the region populated by passively evolving, 
old Es, if they are "moderately attenuated" (i.e., E(B — 
V)star = 0.5) and at 1.3 < z < 2 (Fig. Bl). Conversely, 
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among the model DSGs with a long burst and attenuated 
according to the Calzetti law, no "intruders" are found (Fig. 
B2). Hence the main results of Sect. 4.2 are confirmed here. 

There are no "intruders" with 4 < I c — K < 5 in Fig. Bl, 
as instead in Fig. 8, but this is easy to understand. The orig- 
inal values of the Balmer line-emission reddening E(B — V) 
found by Calzetti et al. (1994) for their sample range be- 
tween and 0.93, so that E(B — V) s tar ranges between and 
0.41 for this sample, iiE{B-V) sta r = (0.44±0.03) E(B-V) 
(Calzetti 1997). Hence the observed values of E(B — V) s tar 
are lower than those considered by PM00 and here. The 
mean of these observed values is consistent with the value 
of E(B - V)star predicted by the WG00 SHELL configura- 
tion with SMC-type dust distributed in a two-phase clumpy 
medium and ry = 1.5 (Fig. 3), which reproduces rather 
well the Calzetti law (Fig. 2) . In Fig. 8 the "intruders" with 
4 < I c — K < 5 have low/intermediate values of ry and, thus, 
E(B — V) s tar ^ 0.3, consistent with the values obtained for 
the Calzetti sample. 

Values of E(B — V) s tar 0.5 (e.g. in Cimatti et al. 
2002a) require either an enormous amount of dust, in the 
two-phase clumpy case, or low amounts of dust and a ho- 
mogenous dusty medium, for the SHELL configuration (Fig. 
3). In both cases the wavelength dependence of the model 
attenuation function is quite different from the Calzetti law 
(Fig. 2) . We conclude that the scaling of the functional form 
giving the wavelength dependence of the Calzetti law with 
arbitrary values of E(B — V) s tar ^ 0.5, under the assump- 
tion that reddening is a measure of dustiness, is not evinced 
from the observations and may be questionable on a theo- 
retical ground (cf. Sect. 2.2.2), as concluded also by WG00. 



This paper has been typeset from a TgX/ P/IpjX file prepared 
by the author. 
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Figure Bl. The distribution of Zq DSGs with a short burst 
and DPSGs in the I c — K vs. J — K colour-colour plane under the 
assumption that their attenuation function follows the Calzetti 
law. Lines are the same as in Fig. 8. Models of selected ages 
are reproduced with different symbols according to their redshift. 
DSGs/DPSGs with a stellar reddening E(B - V) s tar = 0.5 have 
bluer colours than those with E(B — V) s tar = 0.8. DPSGs at- 
tenuated with a Calzetti law "intrude" the region populated by 
passively evolving, old Es, consistent with the result of Sect. 4.2.1. 
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Figure B2. The same as in Fig. Bl but for Zq DSGs with a 
long burst. Models of selected ages are reproduced with different 
symbols according to their redshift. In this case, no "intruders" 
are found, consistent with the result of Sect. 4.2.2 and PM00. 



